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Abstract: We formulate quantum mechanics in spacetimes with real-order fractional 
geometry and more general factorizable measures. In spacetimes where coordinates and 
momenta span the whole real line, Heisenberg's principle is proven and the wave-functions 
minimizing the uncertainty are found. In spite of the fact that ordinary time and spa- 
tial translations are broken and the dynamics is not unitary, the theory is in one-to-one 
correspondence with a unitary one, thus allowing us to employ standard tools of analy- 
sis. These features are illustrated in the examples of the free particle and the harmonic 
oscillator. While fractional (and the more general anomalous-spacetime) free models are 
formally indistinguishable from ordinary ones at the classical level, at the quantum level 
they differ both in the Hilbert space and for a topological term fixing the classical action 
in the path integral formulation. Thus, all non-unitarity in fractional quantum dynamics 
is encoded in a contribution depending only on the initial and final state. 
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1 Introduction 

In the attempt to understand the relation between the ultraviolet (UV) finiteness of dif- 
ferent theories of quantum gravity and their geometric properties, a framework has been 
formulated where the geometry of spacetime displays most of the characteristics of multi- 
fractals, with an anomalous and scale-dependent dimensionality. A traditional perturbative 
field theory lives in a background geometry described by fractional calculus [1-7]. Renor- 
malizability is achieved, at least at the level of power counting, thanks to the change of 
correlation functions with the scale, from a four-dimensional infrared down to small scales 
where spacetime is effectively two-dimensional. This is a general property of models in 
spacetimes with Lebesgue-Stieltjes measures [8-10], which are particular realizations of 



- 1 - 



dimensional flow in quantum gravity [11-13]. Contrary to general Lebesgue-Stieltjes mea- 
sures, the coordinate dependence of fractional measures is factorized [2, 3], which crucially 
allows one to define momentum space and a 'Fourier' transform [5]. Also, fractional mea- 
sures make it possible to realize multi-scale geometries in a quantitative, controlled way 
[6, 7]. 

Depending on its interpretation and applications, multi-fractional theory constitutes 
either a fundamental or effective framework. In the second case, it aims to describe the ge- 
ometry of other models through the tools of multi-fractal geometry, transport and stochas- 
tic theory. These tools are generic enough to allow for a 'portable' description of the physics 
across the models. For instance, the effective cyclic-invariant measures in classical coordi- 
nates of K-Minkowski and other non-commutative spacetimes with non-linear algebras are 
of fractional type and link non-commutative and fractal geometries [4]. Also, the fractal 
properties of the renormalization group flow in asymptotic safety [14] are reproduced via 
fractional multi-scale diffusion equations [15]. 

Taken as a fundamental description of spacetime, on the other hand, fractional geom- 
etry is promising but a number of elementary questions still need to be addressed. A very 
basic one is, in fact, whether quantum mechanics exists in such geometries and, if it does, 
how it is modified with respect to the ordinary case. The presence of a non-trivial measure 
breaking all Poincare symmetries and with a singularity at the coordinate axes may in fact 
carry unforeseen consequences for quantization (first as well as second) and the structure 
of Hilbert spaces and operators. 

The purpose of this paper is to answer the above question in the affirmative and show 
that quantum mechanics on fractional spacetimes is well defined. Thanks to a simple field 
redefinition, the case of a free non-relativistic particle can be formally recast as the one in 
standard geometry, the only difference at the canonical level being in the profiles of the 
wave- functions. As a consequence, even the free theory is non-trivial. It is also very easy 
to work out the quantum fractional harmonic oscillator in all its details. For interacting 
systems (non-linear equations of motion), the underlying genuine fractal/anomalous na- 
ture of space and time is further stressed by the fact that the measure weight cannot be 
reabsorbed completely into the position variables. 

Although ordinary time and space translations are lost, there survives a modified notion 
of them which highlights the existence of a very convenient mapping of the dynamics into 
a unitary one. Even if the dynamics is non-unitary due to the typical dissipative nature 
of fractional models, we shall still be able to introduce self-adjoint operators, manipulate 
probabilities and define scalar products and temporal evolution in almost the same fashion 
as standard quantum mechanics. This is possible for the following reason. In conventional 
dissipative systems, signals decay in time and space and states cannot evolve unitarily. That 
means that there do not exist observables associated with a self-adjoint operator Pq ~ d x 
generating spatial translations and with a Hamiltonian Hq generating time translations. 
The problem is circumvented in the fractional scenario because both space and time are 
modified in such a way that one can construct momentum and Hamiltonian operators 
P ^ Pq and H ^ Hq which are not self-adjoint with respect to the standard scalar product, 
but that are self-adjoint with respect to the natural scalar product on fractional spaces. 



- 2 - 



Thus, one obtains a notion of conserved scalar product and eigenfunctions of stationary 
states do not decay. To achieve this, it is sufficient to notice a one-to-one correspondence 
between operators Aq of the ordinary theory and operators A of the one with factorizable 
(in particular, fractional) measure. In position space, 

i <— >■ A:= }— A ^~(x), (1.1) 

\/Vi{x) 

where Vi{x) is the spatial part of the non-trivial measure weight. For the evolution operator 
U from time t' to t, there is also an extra time dependence in the form \/vo(t')/vo(t), 
where vq is the time-dependent part of the measure weight. Although the theory is non- 
unitary, the iS-matrix limit of the fractional propagator is unitary due to trivialization of 
the term \/ vo(t') / vo{t) in the double limit t' — > — oo, t — > +00, provided vq depends on 
\t\. Furthermore, it is shown that the classical action S c \ evaluated on classical trajectories 
is the same in ordinary and fractional theory. Nevertheless, S c \ enters the path integral 
formulation in a rather interesting way, i.e., augmented by a topological term distinguishing 
the ordinary from the fractional path integral; this term is, again, nothing but a rewriting 
of the ratio y/vo(t')/vo(t). All these features will be discussed in due course. 

After briefly recalling the main properties of fractional geometries and extending the 
discussion to even more general measures (section 1.1), in section 2 we motivate bilateral 
fractional spaces (introduced in [5]) by the interpretational difficulties arising in quantum 
mechanics if both coordinates and momenta are positive definite. We define a pair of self- 
adjoint operators (position and momentum) for quantum mechanics and prove Heisenberg's 
uncertainty principle. The latter, which is formally the same as in the integer case, is 
minimized by a weighted version of Gaussian functions. The Schrodinger equation and the 
operator of time evolution are then constructed. The classical and quantum free particle 
and the harmonic oscillator are studied in section 3, while the path integral formulation 
on fractional spaces is developed in section 4. Conclusions are in the last section. 

Before starting, we comment on other works falling under the unrelated but almost 
homonymous topic of 'fractional quantum mechanics'. It is known that the fine-scale struc- 
ture of a quantum mechanical particle path is very irregular and described by a nowhere- 
differentiable curve akin to a Brownian path [16, 17]. As normal diffusion associated with 
a Wiener process is only one among many possibilities realized by Nature, and anoma- 
lous diffusion described by fractional equations is the rule rather than the exception (see 
[7, 19-21] and references therein), it is spontaneous to ask what type of quantum me- 
chanics would emerge from 'anomalous' paths. Anomalous trasport takes place in chaotic 
quantum systems, of which fractional dynamics is a model in certain regimes [22, 23]. 
Breaking of time reversal and dissipation typical of classical fractional mechanics translate 
to a non-conservation of probability in the quantum realm. These facts prompted curiosity 
in formulating first-quantization dynamics with fractional derivatives. The Schrodinger 
equation was generalized to the case where the Hamiltonian H features a spatial Riesz 

1 This prompted an early attempt to formulate a model of quantum mechanics in fractal spacetime [18], 
not based on fractional calculus. 
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Laplacian [24-32] (stemming from a path integral over Levy paths); long-range interac- 
tions in quantum lattices can be modeled by this type of dynamics [33]. Other modifi- 
cations of the Schrodinger equation involve a non-Markovian time evolution of the form 
[(iH)Pc% — H]i/j = 0, where df is the fractional derivative of order < ft < 1 [34-36] (in 
particular, when ft = 1/2 fractional dynamics can be described by a comb model [37, 38]), 
or both time and space fractional derivatives [39, 40]. See [41] for another type of heuristic 
application of fractional operators to quantum mechanics. 

Our case is quite different with respect to both motivation (as stated above) and 
techniques. The framework we present does not describe anomalous quantum systems 
in an ordinary spacetime but ordinary quantum systems in an intrinsically anomalous 
(in particular, fractal) spacetime. At a loose phenomenological level, this interpretation 
might also hold for some of the cited works, but our technical implementation is unique to 
fractional spacetimes as declared in [1-3, 5]. Morevover, our Laplacian is not a fractional 
operator but a plain second-order operator with fractional weights, chosen from the three 
requirements of (i) reduction to two dimensions in the ultraviolet of a quantum field theory 
[1, 3], (ii) of being a quadratic form and (iii) of self-adjointness [5]. 

1.1 Basic facts about fractional spacetimes 

Fractional spacetimes are a continuum with anomalous properties realizing fractal as well 
as non- fractal types of geometry. They have been introduced in [1, 2] having a field theory 
in mind, so they are formulated as a redefinition of the Lebesgue measure of time and 
space. Here we limit to a minimal set of definitions sufficient to the scope of the present 
paper. In D topological dimensions, fj, = 0, 1, . . . , D — 1 (where x° = t is time), the measure 
dP x = dtd®" 1 *. is replaced by 

D-l 

dg a (x) := d D xv a (x) := dtv ao (t) dx l v ai (xi) , 

i=i 

where < an < 1 are D constants dubbed fractional charges and T is the gamma function. 
There exists also a more fundamental definition of g a based on complex fractional charges, 
which leads to a discrete spacetime texture [1, 3] and is closely related to K-Minkowski 
non-commutative geometry [4], but we shall not consider it here. The real-order measure 
(1.2) is said to be isotropic if all a fjL = a are equal. Without loss of generality, we consider 
a spatially isotropic configuration where all the spatial charges are equal, = a, but the 
time-direction charge ao may differ from them. 

We distinguish between unilateral measures where x M > and bilateral measures where 
the support of v a is the whole real line (minus the axes x^ = 0). In the first case, the 
absolute value in (1.2) is pleonastic and the world is limited to the first coordinate orthant. 

Clearly, the Poincare group is completely broken. In general, one could define a non- 
factorizable Lebesgue-Stieltjes measure weight v(t, x) preserving at least some of the sym- 
metries, as in a previous attempt to realize anomalous spacetimes with a rotation-invariant 
measure [8-10]. However, factorizability of the coordinate dependence is crucial for the 
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definition of the momentum transform in D dimensions, which is [5] 



f(k) := / dQ a (x)f(x)c a>a >(k,x), (1.3a) 
Jo 

f(x) = / dg a >(k) f(k)c aa '(k,x) , (1.3b) 



o 



D 

2\~ 1 



for the unilateral world, and 

/+oo 
dg a (x) f(x)e* a al (k,x) , (1.4a) 
-oo 

/ + CX) 
dg a '{k) f(k)e a ^(k,x) , (1.4b) 
-oo 

^'(M) := , = - — -5-, (1.4c) 

V^a'W^aW (27r)~ 

for the bilateral world. This is only one choice among an infinite class of transforms and 
we have allowed for a momentum-space measure different from position space. Both c a ^ a i 
and <£ a;a ' are eigenfunctions of the fractional Laplacian (in Euclidean signature) 

K- : =E P '' V„:=-fl=dJyfiJx)-} , (1.5) 

n yv a {x) 

with eigenvalue — \k\ 2 . For simplicity we shall concentrate on this second-order operator, 
although operators of fractional order are possible [7]. 

Giving up factorizability in a generic weight v(t, x) quickly leads to a technically un- 
tractable model. Here we shall further stress this point in several occasions with novel 
arguments. Moreover, contrary to non-factorizable measures fractional ones favour con- 
trol over all details, including the rigorous determination of the Hausdorff, spectral and 
walk dimension of spacetime [2, 7] and the formulation of a multi-scale geometry in direct 
contact with the definitions in complex systems [3, 7]. 

Regarding the dimensions, the volume of a Z)-ball with radius R scales as V^' ~ R dn , 
where 

dn = a M (1.6) 

is the Hausdorff dimension [2]. For our particular choice of spatial isotropy, c?h = ao + 
(D — l)a. On the other hand, the spectral dimension ds depends on the type of Laplacian 
and diffusion associated with spacetime. In the simplest case of normal diffusion, ds = dn 
and the walk dimension is exactly equal to 2 [7]. 

The physical picture envisages a spacetime which is effectively two-dimensional at small 
scales, so that field theories therein have improved UV behaviour [3], while at large scales 
one must recover a spacetime with Poincare symmetries restored and the ordinary integer 
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dimensionality dn = d$ = D. To realize this picture of 'dimensional flow,' one needs to 
upgrade to a multi-scale geometry and sum over all possible a's, so that the total measure 
weight is of the form ^ n g n v an for some dimensionful couplings g n [3, 7]. This definition is 
based upon the standard one of multi-fractal geometry. For instance, in an isotropic setting 
the minimum set giving rise to multi-scale behaviour is a £ {a*,l}, where a* = 2/D is 
the UV fractional charge. A direct manipulation of multi-fractional measures can become 
rather complicated, also because they are no longer factorizable and it is not clear whether 
it is possible to define a momentum transform across all the scales [5]. Therefore, it is 
often simpler to work at different 'snapshots' of the geometry by considering a measure 
with fixed, constant charges a^. 

Thanks to the resiliency of the formalism, we can generalize the discussion to factor- 
izable measures which are not fractional, by leaving the Ansatz 

D-i 

dtd ' 1 * — ► dg(x) = dtv (t) JJ dx i v^x 1 ) (1.7) 

i=l 

completely unspecified. At one's own convenience, one can replace fractional measures 
with this general form and vice versa, 

v ao O v , v ai *-> Vi . (1.8) 

As an added bonus, we may even regard Vo(t) as a rau/ii-fractional measure vo(t) = 
9n v a„ (t) limited to the time direction. For simplicity we shall continue to call the 
measures (1.7) fractional, but the reader should keep in mind this important generaliza- 
tion. 

2 Quantum mechanics 

Once defined a unitary momentum transform F g on a suitable L 2 space [5] , it is natural to 
ask whether it is possible to formulate quantum mechanics on fractional spacetimes and, 
more generally, on spacetimes with factorizable measure. We shall: 

(i) define a pair of self-adjoint operators Q and P obeying the Heisenberg algebra; 

(ii) show that Q and P are related to each other by the momentum transform F g , which 
is an isomorphism between two equivalent Hilbert spaces (spanned by x and k); 

(hi) check whether Heisenberg's uncertainty principle holds; 

(iv) find a complete basis of the function space as eigenstates of a self-adjoint operator; 

(v) discuss time evolution, Schrodinger equation and the Green function. 

Applications of the formalism to the free particle and the harmonic oscillator will be given 
in section 3. 
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In unilateral worlds, the momentum transform is real, hence wave-functions are real, 
too, and there is no self-adjoint first-order derivative operator. 2 Moreover, a world with only 
positive-definite momenta would be problematic: particles could only run away from the 
origin, without stopping or coming back, and one would be unable to construct bound states 
and to confine particles in a natural way. This is quite different from the ordinary situation 
of a particle with an infinite potential wall at x < 0, because in that case momentum space 
still would be the whole real line. In quantum field theory, one would face similar paradoxes: 
Feynman diagrams would loose their meaning, as all external legs in a diagram would carry 
only incoming momenta, the concept of anti-particle would be missing, and so on. Looking 
back at the derivation of the unilateral transform [5], it is more sensible to interpret both 
x and k as, actually, the absolute values \x\ and of a bilateral world. Therefore, we 
regard the bilateral world as the correct representation of fractional geometries where one 
can attempt to do meaningful physics. 3 

2.1 Hilbert space 

From here on we specialize to one dimension (and still call Vi the spatial measure weight). 
Let L 2 g (R) be the Hilbert space of functions tjj which are square integrable with respect to 
the non-trivial measure. This space is equipped with the inner product 

/+oo r+oo 
dg(x) ipl(x) ip 2 (x) = / dx Vi(x) ipl(x) ij) 2 {x) , ipi, 2 £ L 2 g . (2.1) 
-oo J — oo 

We shall often employ the bra-ket notation (V'ilV'2) and reserve the symbol {ipi\ip2)o := 
/ dxi/jl(x) ip 2 (x) for the scalar product with respect to the ordinary Lebesgue measure. 
The existence and finiteness of the norm \\ip\\ 2 := (ip, ip) requires that 

V>(x) x ~°° x~ v Re(z/) > | , (2.2a) 

t/;(x) x~ p Re(p) < | . (2.2b) 

Next, we define two self-adjoint operators on this space. For the position operator, we shall 
choose 

Q:=x, (2.3) 

with domain Dq = {ip G L 2 : xip G L 2 g }. Looking at eq. (1.5), it is natural to define as 
momentum operator 



P := —ihV x = —ih — = d x 
\/Vi(x) 



V^) "J • (2-4) 

We establish the domain Dp by requiring P to be self-adjoint. First of all, it will be 
obviously necessary that ^Jvi(x) ijj(x) G C 1 (M), in order to perform derivatives. Notice 



2 Every integration by parts entails a minus factor, so that, at best, a first-order derivative could be an 
anti-hermitian operator. The first available self-adjoint operator would then be second-order, which would 
have the chance to be self-adjoint only because two integrations by parts restore the correct sign. 

3 Since k 2 is positive definite, the derivation of the Green function of a scalar field in [3], worked out in 
the unilateral case, is unaffected. The only difference is that, in a bilateral world, the spectrum is twice 
that of the unilateral one. 
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that the existence of the derivative is required for the function y/vlip, but not for tp. 
Secondly, we should demand that the image of P is in L 2 , Pip(x) G L 2 . Finally, we need 
to check the absence of boundary terms when performing integrations by parts. Let us 
choose for now ip a G L 2 g , a = 1, 2, such that G C^R). Then, 

/+oo 
de(x)ip*(x)'D x ip2(x) 
-oo 
■oo 



-zft / da; Wi(ac) -01 (») 5 X [V ^ 2 (a;)] 

i/ — oo 

+oo r+°° 

-ih[vi(x)i/jl(x)ip 2 {x)] +iH dg(x)V x tpl(x)ijj2(x) 



dg(x) [-iKD x i/ji(x)]* ip 2 {x) 

= (Pih, fa) • (2.5) 

The boundary term in the third equality vanishes at infinity because tpi G Note that, 
due to the presence of the (mild) singularity at x = in the fractional measure Vi(x) = 
v a (x), the behaviour of v a (x)tpl(x)ijj2(x) has to be checked also at the origin: it could be 
either singular or discontinuous. However, since y/v^^a G C 1 (M), then ^(a:)^*^)^^) is 
such (as product of two C 1 (IR) functions) and a fortiori continuous in R and, in particular, 
at the origin. In general, the boundary term indeed vanishes, and with the choice of domain 
Dp = {if; G L 2 : yjvi{x)ip{x) G C 1 (R), Pip(x) G L 2 } the operator P is self adjoint. 

2.2 Position and momentum operators 

In the common domain Dq n Dp, the operators Q and P satisfy the Heisenberg algebra 

[Q,P] = ih. (2.6) 

In a fractional space, the generalized eigenfunctions of the position operator Q are the 
fractional Dirac distribution 5 e , 

Q8 g (q,x) = q8 e {q,x) , (2.7) 

where [5] 

S e (x,x'):= /frT 3 ?,. (2.8) 

is consistent with the resolution of the identity provided by the transform (1.4). 

The normalized eigenfunctions of P, with real eigenvalue p = Hk G R, are given by 
eq. (1.4c). We set F g to be an automorphism (a' = a in fractional measures), which 
simplifies the discussion; a generalization is straightforward. Then, 

P<e g (kx) = p<B g (kx) , e g (fcrr) = — ■= . (2-9) 

yJVi{k)Vi{x) V27T 
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Since the functions <s g {kx) are symmetric under the exchange k •<-)■ x, eq. (2.9) implies also 

1 



y/vi{k) 



dk 



a/ Vi(k) <s e (kx) = xe s (kx). (2-10) 



Equations (2.9) and (2.10) establish the isomorphism between the two equivalent Hilbert 
spaces L 2 {x) and L 2 {k) related by the (unitary) momentum transform (1.4), F e : L 2 (x) — > 
L 2 (k). In fact, once functions are decomposed according to (1.4), eqs. (2.9) and (2.10) 
show that the operators P and Q in L 2 (x) can be equivalently realized as multiplicative 
(by k) and derivative (with respect to the variable k) operators in L 2 (k): in L 2 (k) the two 
operators exchange their roles, and also their domains are mapped one into the other. 4 

It is worth to stress that, for any given measure with a singularity at a point x (in this 
case, x = 0), P (hence H below) is uniquely defined because it is not translation invariant. 
This contrasts with ordinary quantum mechanics, where Po ~ d x is translation-invariant. 

2.3 Heisenberg principle 

The Schwarz inequality, which holds in our vector space, 

|(^1,^)| 2 <||^I| 2 ||^H 2 (2-11) 

is sufficient to prove Heisenberg's principle through the Robertson inequality. Define the 
expectation value of an hermitian operator A as (A) := (i/>, Aijj), and the standard deviation 
about its mean value as (AA) 2 := (A 2 ) — (A) 2 = (A))^|| 2 . For any pair of hermitian 

operators A and B, apply the Schwarz inequality to the vectors ipA = (A — (A))ip and 
tpB = (B — {B))ip, obtaining the Robertson inequality II^aPII^bII 2 > [([A B])/(2i)} 2 . 
When A = P and B = Q, this becomes Heisenberg's uncertainty principle, 

(AQ) 2 (AP) 2 >^, (2.12) 

just like in ordinary quantum mechanics. Any dependence from the non-trivial measure 
has been absorbed into the definitions of the position and momentum operators. 

2.4 Hamiltonian 

We can define a self-adjoint operator corresponding to the fractional quantum Hamiltonian 
of a point particle of mass m in a real potential V(x): 

H = ±P + V(Q) = -£s2 + V (X ) _ [VW)-} + V(«) . (2.13) 



4 It should be noted that the request of having self-adjoint differential operators severely restricts the 
form of the momentum operator in the fractional context and, in particular, it rules out naive fractional 
derivatives. If P were defined by a fractional derivative of some order /3, integration by parts would yield 
an infinite number of terms due a complicated Leibniz rule [2, eq. (2.66)]. Even in the 'trivial' a = 1 case, 
one would still have a further problem: (tpi, d /3 ip2) = (d^tpi, 1P2), where d 13 is the right derivative operator. 
Mixed fractional operators of the form ad 13 + a*d^ exist which are self-adjoint [7], but we do not consider 
them here. 
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For every Hamiltonian of this form, one can always consider the following ordered relation: 



H = -==H oy fc(x), (2.14) 
y/Vi{x) 

where Hq = — [h 2 / (2m)]d 2 +V(x) is the corresponding Hamiltonian operator of the ordinary 
non-fractional quantum theory. By construction, H is self-adjoint. 
The eigenvalue problem 

H^ n (x) = E n il> n (x) , (2.15) 

can be formulated as 

H tp n (x) = E n ip n (x) , (2.16) 

where 



f n (x) := \Zvi(x)ip n {x) (2.17) 

are the eigenfunctions in the standard L 2 space. If the latter are orthonormal under 
the L 2 scalar product, so are the fractional eigenfunctions ip n with respect to eq. (2.1), 

2.5 Time evolution and Schrodinger equation 

As usual in an ordinary quantum mechanical framework, we regard time as a parameter 
and consider the time evolution of the system. The use of a time with non-trivial fractional 
charge is strictly connected to some features of fractional spacetimes regarding violations 
of energy conservation. 

In the Schrodinger picture, the natural generalization of the Schrodinger equation for 
a quantum mechanical state \ip) is 

iKDtlW)) =H\i/>(t)) , (2.18) 

where 



V t = d t 



The choice of T>t in eq. (2.18) is dictated by naturalness (the classical system is defined via 
these derivatives, not dt) and guarantees that quantum mechanics retains the dissipative 
character of the classical theory (a Schrodinger equation defined with dt instead of T>t 
would lead to a unitary model). 5 



5 The imaginary-time version of the Schrodinger equation is sometimes regarded as the definition of the 
diffusion equation valid in a given space, where H is interpreted as the Laplacian K, and t — ► a as an abstract 
diffusion time. The diffusion equation in fractional spacetimes typically features the ordinary first-order 
derivative d a , or a fractional derivative in anomalous models, but never the diffusion operator T> a because 
a is an abstract parameter unrelated from the time fractional structure of the geometry [2, 6, 7]. Thus, 
there is a contrast between quantum mechanics and field theory in fractional spacetimes, marked by a 
non-correspondence between Schrodinger and imaginary-time diffusion equation. Since the difference lies 
in geometric weights, this somewhat reminds one of the situation in general relativity, where the Laplacian 
is the Euclideanized Laplace-Beltrami operator and the diffusion equation is manifestly non-covariant. 
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Let us define, firstly, the quantum state 



|X> := VMt) |V) , (2.20) 

and, then, the unitary state: 

\<p) := y^(xj \ X ) = y/vi(x)v Q (t) |^> . (2.21) 

Notice that Vi(x)vo(t) is the total fractional measure in 1 + 1 dimensions and one could have 
obtained it from the start by defining both T>% and H with a generic measure v(t,x). Fac- 
torizability of the coordinate dependence in the measure plays an essential role in keeping 
time and space dependence separated and it is essential in order to define the eigenvalue 
problem consistently. Indeed, eigenfunctions of the Hamiltonian would depend on time in 
the case of a non-factorizable measure v(t, x), as it is clear by looking at eq. (2.17), and 
would lose their fundamental character as stationary states. 

The state \(p) satisfies the usual (unitary) non- fractional Schrodinger equation 

iW t \<p(t)) = H \<p(t)) , (2.22) 

where Hq is defined in eq. (2.14). Then, considering the complex conjugate of eq. (2.22), 
—ihdt {<p{t)\ = (<p{t)\ Hq, and the requirement Hq = Hq, under the ordinary I? scalar prod- 
uct one gets ihd t (ipi(t)\ip 2 (t)) = (^i(t)|^|^ 2 (i)) - (tpi(t)\H \<p 2 (t)) Q = 0, meaning that 
equal-time scalar products of the type {<Pi(t)\tp2(t)} are time independent, as expected. 

By using the relation (2.20), we see that, in a fractional space, time evolution preserves 
L 2 g products of the form 

(Xi(*)l*2(t)> = (V>i(i)k(*)lV>2(i)> = <^i(to)h(to)|^(to)> = (xi(*o)|X2(*o)) • (2.23) 

This result lets us prefigure that time evolution in a fractional space and the associated 
operator are not unitary, unlike in ordinary quantum mechanics. 

2.5.1 Time evolution operator 

Let Ufa — t') = exp[—(i/h)Ho(t — t')] be the unitary operator associated with the time 
evolution of the state \<p). As it is well known, eq. (2.22) can be written as 

ihd t U (t - \<ptf)) = H U (t - t') \<p(t')) , (2.24) 

and, simply by multiplying and dividing by appropriate measure factors, we get 

1 



ihVf 



y/vi{x)v Q {t) 
1 



Mt-t^^fv^vfa')^ 1 )) 



y/Vi(x)V {t) 

which is the analogue of eq. (2.24) in a fractional space: 

ihV t U{t, t') \^{t')} = HU(t, t') \ip(t')) , (2.25) 
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where U(t,t') is the fractional time evolution operator defined as 

U(t,t') := / e-faWy/viWM*) ■ 

Expanding the exponential operator and resumming, we get a form of U in which the 
fractional Hamiltonian H appears explicitly: 



U(t, f ) := -y=i==Uo(t ~ t>) y/vMM*) = \ VJ m e "^ (t "° 

y/Vi(x)vo(t) V V (t) 



(2.26) 



The ordering of the factors is irrelevant. Notice that we can rewrite this expression as 

U(t,t') = e -iH{t-t')-\!t,dT\nv {T) = e -*S*dT[H-$biMT)] . (2.27) 

The simultaneous presence of non-trivial real and imaginary contributions to the evolution 
operator highlights the dissipative nature of the system. The lni>o term will be further 
discussed in section 4.2. 

One can check that U still keeps some important features characterizing a time evolu- 
tion operator: 

• The initial condition lim^jf \ip(t)) = \ip(t')) is compatible with 

lim U(t, t') = 1 . (2.28) 

t— 

• The composition law required by time evolution holds: 

U{t, t') = U(t, t")U(t", t') , t' <t" <t. (2.29) 

Nevertheless, the operator U shows some novelties characterizing time evolution in a 
fractional space: 

• Time-translation invariance is broken and U is not unitary: 



u(t,t')uHt,t f ) = &(t,t')u(t,t>) = ^ + 1 



so that 
Indeed, 

and 



(2.30) 

U^fj- 1 . (2.31) 



U-\t,t')=ei 6 ^\^ = U(t',t) (2.32) 
V u o(* J 



vo[t) 



(2.33) 



When time-translating the system forwards and, then, moving it back to the initial 
instant, one gains a factor distinguishing the two states. This mechanism avoids the 
so-called 'grandfather paradox' and introduces an arrow of time also at the micro- 
scopic level. 
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• U just implements what we found in the last section, that is, l? Q products such as 
(V'll^ol^) are preserved during time evolution: 

(Mt)\v (t)\Mt)) = (Mt')\uHt,t')v (t)u(t,t')\Mt')) 

= (Mt')\vo(t')\Mt')) • (2-34) 
On the other hand, the scalar product (ipilfa) is n °t preserved in time. 
Thus, while the wave-function 93 evolves unitarily the fractional wave-functions tp do not, 

il>(x, t) = \ [ dt' U (t - t'Mx, t') = [ dt' U(t, t')i;(x, . (2.35) 

^Vi(x)V (t) J J 

2.6 Green function 

As we saw above, the evolution of a fractional quantum mechanical system is described 
by eq. (2.18), ihV t ^(x,t) = Hip(x,t), with an initial condition i/j(x,t) = ip(x,t') at t = t' . 
From a mathematical point of view, we are dealing with a Cauchy problem. Therefore, 
also in the fractional case, the general solution can be obtained by means of the Green 
function G(x,t;x',t') and, in particular, by solving the following system: 

(ihV t G(x, t; x',t>) = HG{x, t; x',t') , 
{G(x,t';x',t') = 6 e (x,x') . 

The solution of eq. (2.18) can be represented as 

ip(x,t) = J dx'vi{x')G(x,t;x',t')i;(x',t') . (2.37) 

Following the same route outlined above, we rewrite eq. (2.36) as iHdtGo(x, t; x', t') = 
HqGo(x, t; x', t'), where Go := >/vi(x)vo(t)G and ^0 are the unitary quantities of the 
ordinary quantum theory. The solution of the last equation is well known: 



G (x,t;x',t') = c^2<p n (x)ipn(x')e a 



E n {t-t') 



where c is independent of x and t, and (p n are the eigenfunctions of the hamiltonian Hq. 
The value c = 1 corresponds to the initial condition Gq(x, t'\ x', t') = 8{x — x'). The form 
of G is then 

G(x, t; x', = , , C , = E • (2-38) 



By imposing the initial condition of eq. (2.36), we get the constant c = y / vo(t')/vi(x') and 
we can finally write the fractional Green function as 



G(x,t;x>,t>) = J^^MxWnWe-i*"™ 



(2.39) 



where we used the relation (2.17) and ijj n are the eigenfunctions of the fractional Hamil- 
tonian H. Again, factorizability of the measure is fundamental. Using the completeness 
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relation of the eigenfunctions, it is straightforward to reexpress eq. (2.39) via the time- 
evolution operator: 

G(x, t; x',t') = (x'\U(t, t')\x) . (2.40) 

All the non-unitarity of the theory is encoded in the pre-factor \/vo(t')/vo(t). In section 
4.2 we shall provide a neat reinterpretation of this pre-factor as a topological term in the 
path integral formulation. For the time being, we notice that the S matrix, defined as 

lim (x'\U{t,t')\x) , (2.41) 

t'— CO 

is unitary in fractional theories, provided the limits are performed symmetrically, so that 
\/ v a (t')/v ao (t) — > 1. This is not true in general for theories with measures (1.7), even if vo 
is manifestly positive. An example is vo(t) = exp(at), where a is a constant. A necessary 
and sufficient condition is vq to depend on \t\. 



3 Examples 

In this section we shall develop the two basic quantum mechanical examples: the free 
particle and the harmonic oscillator. 

Formally, a classical mechanics system is insensitive to the choice of spatial measure. 
Its fractional action is 



S 



t-2 



dtvo(t) L . 



(3.1) 



where only vq appears. In order to get the equation of motion, we apply the variational 
principle by finding the stationary points of the action, 



5S = 0, 



(3.2) 



under variations 5q(t) preserving the extrema Sq(ti) = Sqit^j = 0, namely, 



dtvo(t) 



ti 



dL ^ dL 



Here we assumed that the Lagrangian depends only on q and V t q. Since 



(3.3) 



5V t q = V t (q + 5q) - V t (q) 
1 



1 



-.0, 



Vvo(t) 
V t 5q , 



dt 



\fvo{t)5q 



(3.4) 



integrating by parts eq. (3.3) we have 



1 2 



dtvo(t) 



T>1 



dL dL 



dV t q dq 



5q(t) + 



Mt)^ q 6q(t) 



(3.5) 
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The last boundary term vanishes by virtue of 8q(ti) = 5q(t2) = and one is left with the 
fractional generalization of the Euler-Lagrange equation of motion 

n,«_«_0. (3.6) 
oD t q oq 

Moving to Hamiltonian formalism, we define the conjugate momentum 

which allows us to write the Hamiltonian as 

H:=pV t q-L. (3.8) 

Upon postulating the canonical brackets {q,p} = 1 (the others vanish), fractional Hamil- 
ton's equations are 

T> t q:={q,H}, V t p:={p,H}, (3.9) 
that are easily seen to be consistent with (3.6). 

3.1 Free particle 
3.1.1 Classical particle 

The simplest example is that of a free particle. Its classical fractional Lagrangian is defined 
as 

L = -m{V t qf , (3.10) 
leading to the Euler-Lagrange equation of motion 



V\q = ~^=d 2 t \VMf)o\ = • (3.11) 



Notice that by introducing the variable 6 



X (t) = VMt)q(t) , (3.12) 

the fractional equation of motion written in terms of \ would be undistinguishable from 
that of a free particle in a non-fractional space, i.e., 

X = 0. (3.13) 

The classical fractional solutions are then easily written in terms of the non-fractional 
solutions x(t) = Xo + wt, leading to 

q(t) = ^=T=e{xo + wt). (3.14) 



3 This setting is the classical analog of the quantum state \x) introduced in (2.20). 
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Since xo is an arbitrary constant, the first term of this equation highlights a symmetry of 
the Lagrangian (3.10) in a time-dependent change of the position, q(t) — > q(t) + Xo/\/ u o(*) ; 
which reduces, for fractional measures, to an ordinary translation when a = 1. 

The classical fractional action evaluated on the solutions of the equation of motion is 
identical to that of the corresponding non-fractional system x> 

S* c i = jT 2 dtv (t)^m(V t q) 2 = ^mw 2 {t 2 - h) . (3.15) 

In spite of this fact, the interpretation is different: w is not the velocity (which is not 
constant) of the fractional variable q(t). 

The redefinition (3.12) and the corresponding trivialization of the problem obviously 
triggers the following remark. Since this (and also the harmonic oscillator below) is a free 
system, the measure weight can always be reabsorbed into the definition of the coordinates 
and the physics is classically equivalent to the usual one identically. 

Therefore, there is seemingly no point in working out fractional rules (variational prin- 
ciple, time evolution, and so on) which turn out to be unnecessary, and fractional dynamics 
seems altogether trivial. Even for a non- linear system with a potential V(q), one could sim- 
ply consider an ordinary setting with a non-autonomous potential f(t)V(x), for some func- 
tion /. For instance, a power-law potential would yield vo(t)V(q) = vo(t)q n = n ^ 2 {t)x n - 
Are fractional spacetimes just trivial reformulations of non- autonomous systems? 

At the classical level, the answer is mainly interpretational. A non-autonomous poten- 
tial defined in an ordinary space can be considered by all means, but a justification for a 
particular choice should then be provided. So, while in ordinary spacetime a potential of 
the form \t\ l x n with I > would come out of the blue, its origin as an effective description 
of a system with fractional time and the subsequent relation I = (n/2 — 1)(1 — oq) between 
physical parameters are clear. In turn, fractional time and the specific form of the time- 
dependent parts of the system is grounded upon a rigorous implementation of fractal (in 
general, anomalous and multi-scale) geometry via fractional calculus [2, 3, 5, 7], along the 
motivations outlined in the introduction. The multi-scale anomalous/fractal picture and 
its tools cannot be forfeited without loosing logical coherence and the motivational drive. 

However, this perspective is not particularly compelling in the case of free systems, 
where one can erase all measure dependence and work with x coordinates. At the classical 
level, one could simply bypass the problem by adopting a Laplacian of fractional order [7] 
and invoking again the interpretation according to fractal geometry. Unfortunately, this 
would not amount to an actual solution, because second-order Laplacians are desirable in 
a field-theory setting [3]. Quantum mechanics comes to the rescue with a twofold simple 
answer. First, while the dynamics is formally reduced to the standard case due to the 
linearity of the Schrodinger equation, even free systems are not equivalent to non-fractional 
ones because the Hilbert space and the profile of wave-functions drastically change. This 
will be made clear in the example of the harmonic oscillator. Secondly, at the level of the 
path integral, the key difference between ordinary and fractional systems is a topological 
term, a total derivative immaterial in the classical theory. 
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3.1.2 Quantum particle 

The quantum mechanical system can be easily obtained starting from the classical Hamil- 
tonian that, according to (3.8), is H = p 2 /(2m). A self-adjoint Hamiltonian operator H 

-2 
J Q 



in L 2 is promptly obtained by replacing p with the self-adjoint operator P as in (2.4). 



Then the eigenfunctions of the Hamiltonian are just (2.9), with eigenvalue h 2 k 2 /(2m). The 
spectrum is continuous, and the eigenfunctions i^k( x ) = <£ g {kx) are normalized per unit 
volume. 

The Green function of the system can be obtained from (2.39) (with the obvious 
replacement ^ n — > J dkvi(k) due to the continuous spectrum): 



G(x,t;x',t' 



^ [ dkv^^ikxXJkx^e-^ 2 ^ 



vo(f) 1 



v (t)vi(x)vi(x') 2ttH 



(t)vi(x)vi(x') V2nh(t-t') ' 1 • ) 



where, to perform the Gaussian integral, we used the regularization m — )■ m + ie and 
performed eventually the limit e — > 0. The argument of the exponential in the last equality 
of (3.16) is just (i/fr) times the action (3.15), that is ew R ' Scl . Note, in fact, that the 
corresponding non- fractional \ particle would have constant velocity w = {x — x')/(t — t'). 
Later we shall come back to the significance of the pre-factor \/vo(t')/vo(t). 

3.2 Harmonic oscillator 

As another application, we examine the harmonic oscillator. Firstly, we build the classical 
fractional theory and obtain the equation of motion as well as the dynamical variables. 
Secondly, we move to the quantum theory, getting a complete orthonormal basis of L 2 by 
solving the associated eigenvalue problem. 

3.2.1 Classical oscillator 

We define the fractional action of the harmonic oscillator of mass m and frequency u as 



rt2 rti 
S = dtv (t) L= dtv (t) 
Jti Jti 



-m (p t q) - -muj q 



(3.17) 



so that one gets the Euler-Lagrange equation 

V 2 q + co 2 q = 0, (3.18) 

that is 



,+!5 4+ if--; ! l+« , ]« = o. (3.19) 

v 2 \v 2vfi 



In the fractional case vo(t) = \t\ a ° 1 /r(«o); we recognize an anti- friction and a potential 
term: 



1 - q . 



(1 - ao)(3 - a ) 2 

77o r u 

At 2 



0. (3.20) 
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For olq = 1, we recover the standard non- fractional equation. 

Moving to Hamiltonian formalism, we define the conjugate momentum 

dL 

p: =m- q =mVtq > (3 - 21) 

which allows us to write down the Hamiltonian as 

H := pV t q — L = —p 2 + -muj 2 q 2 . (3.22) 
2m 2 

Fractional Hamilton's equations turn to be 

V t q:={q,H} = , V t p := {p, H} = -moj 2 q , (3.23) 
m 

where the definition of time evolution is consistent with the result (3.18). Invariance of the 
dynamics under the exchange of p and q variables is evident in the case oj = 1/m. 

Like in the free-particle case, an equivalent way to get the same results is by defining a 
new coordinate variable % := y/voQ an d a new Lagrangian L$ := vqL. The dynamics is then 
the usual (non-fractional) one. The equation of motion one gets is indeed x + ^X = 0) 
which is nothing but eq. (3.18). The x solutions are x(t) = A cos cot + Bsinojt, from 
which the solution q(t) trivially follows. By choosing boundary conditions xi = x(ti) anci 
X2 = xih), the x(t) solution reads 

X(t) = . r 7- - { X 2 sm[u(t - h)} - xi Bm[u(t - t 2 )]} , (3.24) 
sm[o;(t2 — t\)\ 

and the classical action evaluated on the solution (3.24) reads 

Scl = 9 . f ™? 7 y\ {(Xi + Xl) cos[^fe - h)\ - 2xiX2} ■ (3.25) 

Just like in the free particle case, the classical action of the x system (3.25) is undistinguish- 
able from the classical action So we would have obtained substituting the corresponding 
fractional solution q(t) = x{t) / \J vo(t) in the fractional action (3.17). 

3.2.2 Quantum oscillator 

A complete orthonormal basis of L 2 can be obtained by solving the eigenvalue problem 
of the fractional harmonic oscillator. Set h = 1 and also m = 1 = u). The following 
Hamiltonian is manifestly self-adjoint as sum of products of self-adjoint operators: 

H = \ P 2 + \Q 2 . (3.26) 



The eigenvalue problem 

1 



- K 



y/vi(x)i/j n (x) +^x 2 if) n (x) = E n ip n (x) (3.27) 
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is trivially solved once recognized that eq. (3.27) is just the usual harmonic oscillator 
eigenvalue equation when written in terms of the function <p n (x) = y/ Vi(x)tp n (x) , with 
normalized solutions and eigenvalues 

-| d%(p n (x) + \x 2 (p n (x) = E n ip n (x) , 

/ x e~~ H„(x) „ 1 „ „ 

<p n (x) = —r==4 > E n = n+-, (3.28) 

where H n are Hermite polynomials: 

H n {x) := (-l) n e x2 d™e- x2 , n G N . (3.29) 

The eigenvalues E n and eigenfunctions ip n can be also obtained via the method of ladder 
operators; thanks to the self-adjointness of Q and P, in fact, it is possible to define consis- 
tently a set of creation, annihilation and number operators. The proof is identical to the 
standard case and we omit it. Since {(p n (x) \ n £ N} is a complete orthonormal set in L 2 , 
then 

wo = = ;7^» , (3.30) 

is a complete orthonormal set in L^(x), and the Vvi's are the eigenfunctions of the Hamilto- 
nian (3.26) with eigenvalues E n = n + 1/2. The Hamiltonian (3.26) is left invariant under 
the exchange P Q, so the corresponding equation in momentum space is identical, and 
its solutions are therefore the same functions (3.30) with x replaced by k: this means that 
the orthonormal set (3.30) is left invariant by the momentum transform F g . 

The wave-functions minimizing the uncertainty principle (2.12) are weighted Gaus- 
sians: 

(x-x Q ) 2 



i>{x)= === , ||^l| 2 = l. (3-31) 
■x/vi(x)^ir 

Since (Q) = x , {Q 2 ) =x 2 Q + 1/2, (P) = 0, and (P 2 ) = 1/2, we obtain (AQ) 2 (AP) 2 = 1/4. 
Note that for xq = it gives just the lowest eigenfunction iPq{x). In fractional spaces, 
the factor [u j(x)] -1 / 2 = [v a {x)}~ 1 ^ 2 oc |x|^ 1_Q ^ 2 guarantees that the probability to find the 
particle at the origin, where the measure is integrably divergent, is zero. 

It is important to stress that the physics is not translation invariant because of the non- 
trivial measure and the shifted wave-function (3.31) ip(x + &*), where is some constant, 
no longer minimizes the uncertainty principle. This is not a surprise: translation invariance 
is broken in fractional theories. In spite of this fact, all the functions (3.31) saturate the 
Heisenberg bound, for any value of xq. This means that, although translation is no longer 
a symmetry, some similar symmetry survives. Indeed, it is the symmetry generated by the 
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Figure 1. The fractional wave-functions (3.31) minimizing the uncertainty principle, peaked at 
xq = (thick curve, left panel) and near xq = 1 (thick curve, right panel), for a = 1/2. The dashed 
curves are the Gaussian in the ordinary case, a = 1. 



P operator. Let us consider its action on the basis functions ip n = cp n (x)/ \/vi(x): 

(-x ) m P m <p n (x) 



ex.p{-x P)tp n (x) = ^2 



m=0 

(-x ) m 



0': 



/—^ exp[-x d x ](p n (x) 
V v i( x ) 

if n {x - Xq) 



(3.32) 



'P-translations' work in a very convenient way: they first remove the translation non- 
invariant (and non-analytic) factor [^(x)]" 1 / 2 , then act as a translation, and then put 
back the removed non-analytic factor. In other words, assuming that Uj always depends on 
\x\, P-translations shift only the analytic part of the wave-functions. Thus, the system is 
not invariant under spatial translations but it is invariant under P-translations. The eigen- 
functions of the fractional harmonic oscillator equation are t(? n (x; xq) = <p n (x — Xq ) j \Jv{ (x) , 
for any xq, so that the Gaussians (3.31) saturating the Heisenberg bound are nothing but 
the most general lowest-level eigenfunctions tpo(x; xo) (including P-translation degeneracy). 

In the ordinary case, translational symmetry trivially corresponds to different choices 
of the origin of the coordinate system. In the fractional case, different choices of xq give 
different functions. In figure 1 the wave-function (3.31) is depicted for m = v a and the 
values xq = and xq = 1; it vanishes at the origin, where it is not analytic (cusp). Contrary 
to the ordinary case, there are two local peaks, the one near xq being the higher. When 
xq = 0, both peaks have the same height and are symmetric with respect to the origin. 
These are two of the lowest eigenfunctions (parametrized by xq) of the fractional harmonic 
oscillator, and both saturate the Heisenberg principle. 

The presence of cusps is a clear evidence of non-analyticity of the wave- functions, which 
are not differentiable at the origin. However, ordinary derivatives are not the most natural 



differential operators in this fractional context, as they never correspond to hermitian 
operators. Only the combination defining P is hermitian and, regardless the presence of 
cusps, all the harmonic-oscillator wave-functions are C°° with respect to the P-derivatives. 
This property is somewhat reminiscent of the fact that most fractal curves and other 
functions common in chaos theory are nowhere differentiable in the ordinary sense, but 
they can be differentiable in the sense of fractional calculus (e.g., [42, 43]). 

Once given the complete set of normalized eigenfunctions, we can compute the Green 
function using eq. (2.39). To this purpose, one should first recall that, for any \z\ < 1, [45, 
eq. 18.18.28] 

oo .. 

V M*)My)z n = =e [4xy Z -(x*W)(iW))/[2(i- z >)}_ (3 _ 33) 

n=0 ^TTVi(x)Vi(y)(l - Z 2 ) 

Then, applying the above formula with z = e — )/(fr+* e ) (i n such a way that \z\ < 1) and 
eventually performing the limit e — > we get 



G(x ,._/ t >) _ / "q(*0 / 1 e i2iI^[(* 2 +*' 2 Wt-t')-2^'] 

{,h ' ' Y VQ{t)vi{x)vi{x>) ]j 2mhsm(t-f) 

(3.34) 

As in the free particle case (usual and fractional), the argument of the exponential is the 
classical action evaluated on the classical trajectory. Again, however, it should be stressed 
that the fractional action evaluated on the classical trajectory q(t) is undistinguishable 
form the classical non-fractional action of the corresponding x-labelled system. This means 
that, as it is, the exponent in (3.34) is unable to discriminate between fractional and non- 
fractional systems. All the difference is in the term y/vo(t') /[vo(t)vi {x)vi{x')\. 

4 Path integral 

The path integral can be formulated in a fractional context. The Green function previously 
obtained shall be recovered as a path integral. This formalism will permit to clarify the 
role of the non-unitary pre-factor of the fractional theory, providing an interpretation of 
the latter clS el topological term. 

4.1 Quantum mechanics from path integral 

As is well known, the path integral picture of quantum mechanics relies on two postulates. 
It is worth repeating them to understand how these should be intended in a fractional 
world. 

1. The Green function G(x,t;x' ,t') of the system is a sum of phases, 

G(x,t;x',t') = Cj2 e^^'l (4.1) 

paths 

The righ-hand side of this equation is entirely classical. It is a sum over all possible 
classical trajectories joining the point x' at the time t' with the point x at the time t. 
The argument of the exponential is the classical action evaluated at that particular 
trajectory. All the trajectories have the same weight. 
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2. C is a constant to be determined in such a way that G behaves as a probability 
density: for any intermediate i € one must have 

G(x,t;x',t') = J2G(x,t;x,i)G(x,t;x',t'). (4.2) 

Let us begin from (4.2). If the space has anomalous measure g(x), the sum over the 
intermediate positions x should be interpreted as an integral with measure dg(x) = dx Vi(x). 
But then, the identity Green function (corresponding to no propagation) should be realized 
by a fractional delta, eq. (2.8). In turn, the identity Green function is attained for t = t' 
and then i = t' in eq. (4.2). Consequently, this equation can be interpreted as a resolution 
of the identity when evaluated at t = t = t/. Comparing this with (1.4), whose integral over 
dp provides a resolution of the identity, one can conclude that all the (x,x') dependence in 
C is of the form 

C = C (*'*) , C(t, t) = const. (4.3) 
^/vi(x)vi(x') 

where C(t', t) is a function that can depend on time both explicitly and implicitly, through 
vo(t). The latter dependence can be easily recovered by simple scaling arguments. From 
now on we shall specialize to the case of quadratic Lagrangians, so that all the vo(t) factors 
can be reabsorbed with the redefinition = \/vo(t) x(t), while the whole x dependence 
in the action is through x an d its time derivatives. For instance, in the free-particle case 

^mj drv (r) [V t x(t)} 2 = \m j dr [ X (r)} 2 . (4.4) 

Consequently, in the small t — t' limit the resolution of the identity can only provide a 5(x)i 
which scales as S(x)/y/vo(t). Hence the implicit t dependence of C(t',t) is 1/y vo(t). In 
turn, this fixes the vo(t') dependence as y/vo(t'), in such a way the dependence cancels 
at t = t'. Summarizing, initial conditions and scaling arguments fix the following form of 
C: 

C = K \S f,WW^ < 4 ' 5 > 

V V (t)Vi(x)Vi(x') 

where K is constant both in space and time. 

Although most of the considerations will hold for any quadratic Lagrangian, we shall 
continue for simplicity with the free-particle case. The sum over all the possible paths in 
(4.1) can be conveniently performed by slicing the time interval t — t' into iV infinitesimal 
parts, 

t = t N > t N -i > t N - 2 > • ■ ■ > t 2 > h > t = t' , t n -t n -i = e, Ne = t-t', (4.6) 

and then enacting the standard discretization of variables x n = x(t n ). Notice that this 
is the ordinary splitting except for the requirement that none of the points coincide with 
the singularity of the measure; in the fractional case, t n ^ 0. While this constraint on 
the internal points is immaterial due to the eventual analytic continuation leading to the 
final result, it still holds for the initial and final points t and t'. The singularities in 
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the measure must not correspond to either boundary point. Apart from this, there is no 
other complication and we do not even need to redefine time as in other path integrals 
featuring integral singularities, like the hydrogen-atom system [44]. We shall evaluate first 
the discretized Green function Gn(x, t; x' , t'), performing eventually the N — > oo limit. The 
constant K in (4.5) is also defined cis cL limit, K — lirriTv—^Qc ktv, kjv being the corresponding 
constant in the discretized Green function which, by virtue of (4.2), can only be a pure 
power, Kjy = (k) n ■ 

Within the n-th slice, the action attains the value 
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im 
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(4.7) 



Iterating N — 1 times eq. (4.2), we obtain 



Gn(x, t; x' , t') 



v (t N )vi(x N )vi(x ) 
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j dx\ ■ ■ ■ dxN-ie Tl Sn 
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Vn V v o(tN)vi(x N )vi(x ) 
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where we set xo = v^ol^o) = x' an d Xn = ^/uo(ijv) %N = X- Note that all the V{ factors 
cancelled pairwise, except the first and the last. The integration measure simplifies, and 
all the integrals are Gaussians. Positivity of the measure is instrumental to perform all 
simplifications without picking up phases. In the last equality we performed the change of 
variables \n = \/vo(t n ) x n , with Jacobian 
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N-l 
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i 



V N n=1 VMtn) 



(4.9) 



Here Vat is the regularized (discretized) determinant of the operator defined by the kernel 
V(t, t) = vo(t)6(t — r) in the interval (to,tjy). Its N — > oo limit is either 1 or an ill-defined 
constant (zero or infinity) that will be reabsorbed in the constant k n prior to the N — > oo 
limit. The remaining Gaussian integrals are straightforward, leading to 

G N (x,t;x' ,t') = G N (x N ,t N ;x ,t ) 



2 2ireh \ 



vm 



vo(k 



imV l J[ N ) V 2-KhNe y v {t N )vi(x N )vi(x ) 



e 2fijv e . (4.10) 



In performing the limit N — > oo in eq. (4.10), one should recall that Ne = t — t', so 
that the argument of the exponent is just the action of the classical non-fractional particle 
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evaluated on the classical solution (a particle with constant velocity w = (x ~ x')/(t ~ 0) • 
The only delicate term is the first one. Its large N limit is undefined, unless 



N / ■ \ N/2 

( im ^ . (4.11) 
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For such a value, the first factor gives 1 for any N, and the Green function of the free 
particle exactly coincides with eq. (3.16), 



G(x,t;x',t'):= lim G N (x,t;x> ,t') = J V f] J ™ e^W^T, (4.12) 



where, we recall, x = \/ v o(t) x an d x' = yVoW) x ' ■ ^ n similar fashion, it can be shown 
that also the Green function (3.34) of the fractional harmonic oscillator can be obtained 
along the same lines. 

4.2 Non-unitarity and topological considerations 

Whatever the method chosen to evaluate the Green function (either via eq. (2.39) or 
through the path integral), we encountered the following general property: all the non- 
unitarity of the theory is encoded in the pre-factor \/vo(t')/vo(t) appearing in eqs. (3.16) 
(or (4.12)) and (3.34) and in the general expression (2.39) (the Vi terms, instead, provide 
the correct normalization factors for the L 2 g space) . Were it not for such a term, the theory 
would be unitary and indistinguishable from the corresponding ordinary one. As we are 
going to see now, understanding the origin of this pre-factor also clarifies the main dif- 
ference between ordinary and fractional theory at the quantum level, and the reason why 
classical free theories are equivalent. 

The non-unitary term can be rewritten as an additional exponential of a time integral 
(b 7^ is an arbitrary constant) 



V V (t) 

= g- It' drd r ln l b \/ Vo(t)] 

= e -l drV T (b) _ ( 413 ) 

In the non- fractional case, the Green function of a free particle is the exponential of (i/h)So, 
So being the classical action evaluated on the classical solution. In the fractional theory, the 
classical action S c \ evaluated on the classical solution is identical to the free case, 5 c i = 5o 
(see eqs. (3.16) and (3.34)), so the same classical action cannot serve as exponent of the 
quantum Green function also in the fractional context. The crucial and only discriminator 
between fractional and ordinary free theory is the total derivative (4.13). In fact, the 
non-unitary pre-factor needed to obtain the correct Green function can be completely 
reabsorbed in a redefinition of the action: 



1 



S c i — ► S qu = S cl - - I drV T (b) . (4.14) 
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The difference between S c \ and S qu is a boundary term and the two actions are classically 
equivalent (their difference is an integral of a total derivative that does not contribute to 
the equations of motion). Quantum mechanically, they are not: both of them describe free 
theories, but the former in a standard spacetime, the latter in a fractional one. Therefore, 
all the effect of a fractional spacetime can be encoded in a purely topological term. 

To summarize, the traditional recipe G ~ e ( i A) ,s b f or the path integral is incomplete 
inasmuch as two systems having the same classical action on classical trajectories may 
be still physically inequivalent quantum mechanically, and they can be distinguished by a 
topological term. This term strongly reminds us of a 1-cocycle [46], although its interpre- 
tation is different. While 1-cocycles arise whenever the action possesses a symmetry not 
enjoyed by the Lagrangian and they enter as projective representations of such symme- 
try, in this case the topological contribution tells apart ordinary from fractional quantum 
systems. 

Incidentally, this suggests an 'infinitesimal' version of the coordinate transformation 
Q = x/V^o mapping ordinary to fractional classical free systems (and vice versa, but here 
we analyze only one way). Consider an ordinary quadratic Lagrangian and a coordinate 
transformation \ ~ * X + $X- On shell, the change in the Lagrangian is 

9L . dL r _ 
OX OOtX 

dt(p x Sx) + A(l)Px 5 X 

V t (p x S X ). (4.15) 



(3.13) 



On the other hand, eq. (4.14) implies that to generate the topological term and obtain the 
fractional classical action one should impose 5L = —T>t(l) which, compared with (4.15), 
yields 

Sx = ~— ■ (4-16) 

Px 

This is the infinitesimal transformation mapping ordinary into fractional systems. 
4.3 Path integral from quantum mechanics 

The purpose of this section is complementary to the previous ones: we wish to recover the 
path integral postulate (4.1) entirely from a quantum mechanical setting. We fix m = 1. 

Let us follow the standard way by considering a time interval and subdividing it into iV 
small parts as in eq. (4.6). Recalling eq. (2.40), we introduce N — 1 completeness relations 



dxvi(x) \x) (x\ = 1 . (4.17) 

Then, the time evolution operator must be factorized as 

U(t,t') = —L=e-iK-i^...e-i^y^(F), (4.18) 
y/vo[t) 
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and it is clear that the following formula holds: 

. N-l 



I , t ,s nN-1 N-l 

G N (x,t;x',t') = a / ] dx 8 Vi(x s ) T\ G n+ i n (x n+ i,x n ,t n+ i,t n ), (4.19) 



where 
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Consider a Hamiltonian operator of the form H = P 2 /2 + V(Q). Then, 



{x n+1 \H\x n ) = (x n+1 \P 2 + V{Q) 
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where we used eq. (2.9) for the eigenfunctions of the fractional momentum operator P. 
Thus, we have 
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27rh V 'vi(x n+1 )vi(x n ) 
The various spatial measures in eq. (4.19) simplify except for those corresponding to 
(x,t) and (x',t'), and we have 
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Integrations over momenta are Gaussian, providing 
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where, in the last equation, we used eq. (4.11) and changed variables x n = \/vo(t n ) q n . 
The Jacobian of the change of variables is precisely that in (4.11), and we used the fact 
that the potential is quadratic to scale a vo factor out of V. 

We now perform the large N limit to get the Green function G(x, t; x' , t'). In this limit, 
the first factor gives (4.5), the differentials Y\ n dq n provide the definition of the functional 
measure [dq], and the action in the exponential is just the fractional action, 

\m{V t qf - V(q) }, 

(4.25) 

in full agreement with section 4.1. For non-quadratic Lagrangians the problem is subtler 
but, in those cases, the path integral is typically unsolvable. 

Notice that the last line of eq. (4.23) defines the relation between Hamiltonian and 
Lagrangian. The discrete expression eLo,™ = Pn+i{x n +i — x n ) — eH n , H n := p^ +1 /(2m) — 
V(x n ), translates into the Legendre transform Lq = px — H, where Lq is the ordinary 
Lagrangian. Changing variables before integrating over the momenta, this relation becomes 
L = pT>tQ — H for the fractional theory, in agreement with the Ansatz (3.8). 

5 Discussion 

Spaces with non-trivial but factorizable measures provide an interesting arena where to 
explore avenues towards ultraviolet-complete models of Nature. All questions asked in 
standard scenarios of space, time and matter should be addressed also in these anomalous- 
spacetime models. As an application of the existence of unitary invertible transforms [5] 
between configuration and momentum fractional spaces, we have constructed a well-defined 
quantum mechanics on such spaces, proving Heisenberg's principle, formulating the path 
integral and considering the standard examples of the free particle and the harmonic os- 
cillator. In the latter case, the energy levels of the Hamiltonian are the usual ones but 
the eigenfunctions of the problem are augmented by a factor dependent on the Hausdorff 
dimension of the space. These wave-functions have several unusual properties. They are 
i.e., non-analytic in the origin but everywhere differentiable, infinitely many times, 
with respect to the first-order fractional momentum operator P. This operator gener- 
ates a modification of ordinary translations which do not leave invariant the shape of the 
wave- functions, yet it transforms one into another with the same Heisenberg uncertainty. 
Similarly, the time evolution operator realizes a fractional modification of ordinary time 
translations. In general, the theory is non-unitary due, however, to a topological effect de- 
pendent on the initial and final state. This effect, which tells apart fractional from standard 
quantum dynamics, clearly arises from a topological term in the path-integral formulation. 

Before concluding, we examine an issue related to the fractal character of quantum 
mechanical paths. As we recalled in the introduction, the path of a free quantum particle is 
continuous but nowhere differentiable [16]. This property, together with their self-affinity 
(in general, these curves are not self-similar), counts quantum paths as random- fractal 
curves. At scales larger than the Compton wave-length of the particle but smaller than the 
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de Broglie wave-length, in D = 1 the graph (i.e., the plot of q(t) versus t) of a quantum- 
mechanical non-relativistic particle has Hausdorff dimension dn = 3/2, while the path 
(i.e., the 'trace' left by the particle in configuration space q) has dn = 1. In D > 2, 
both have dn = 2 [17, 47-53]. At smaller scales du = 1, as one can understand from 
the analogy between quantum paths and Brownian motion [48, 50, 54]. The relativistic 
particle was studied as well, and dn = 1 [55]. Given the intrinsic 'fractal' nature of the 
classical spacetime wherein the particle moves in fractional theories, one could wonder how 
the above results on the Hausdorff dimension of graphs and paths are altered. The answer 
is that they are not altered at all. One can show this by an explicit calculation following, 
e.g., the one in [17] based on wave-functions and the time evolution operator. One starts 
from the operational definition of fractal dimension (which, in most cases, coincides with 
the Hausdorff dimension dn)- Since fractional spacetimes are continuous, this 'fractal 
dimension' is indeed dn and we can simply refer to it as the Hausdorff dimension. Let Ia 
be the ordinary length of a curve measured with a probe with some given finite resolution 
Ax 0. The Hausdorff length Zh of the curve is defined as In — h m Az->o Za(Ax) H_1 . For 
a smooth curve, dn = 1 and the Hausdorff and ordinary lengths coincide, while for fractal 
curves only In is a meaningful geometric indicator. In fact, as the resolution Ax decreases 
I a tends to infinity (because one needs, say, more and more balls to cover the curve); the 
Hausdorff dimension is always larger than the topological dimension of the fractal set (in 
this case, l). 7 Requiring that In be resolution independent yields the correct value of the 
Hausdorff dimension; the example of the von Koch curve can be found, e.g., in [53]. 

For a quantum particle, I a is replaced by the average distance (I a) covered by the 
particle in a time interval At = t — t' . The average is calculated on a state associated 
with a sequence of position measurements localizing the particle in a region of size Ax. 
The only difference, immaterial for the final result, with respect to the standard case is 
in the non-unitary pre-factor vo(t')/vo(t) in the expression for {Ia)- Up to this factor, 
clearly inherited from the fractional operator U (appearing twice in the expectation value), 
the scaling of the average distance is the standard one, (I a) °c At /Ax. Replacing this 
into In and demanding the Hausdorff length to be resolution-independent yields dn = 2. 
This result was expected because the quantum dynamics of the free particle is the same in 
ordinary and fractional spacetimes except for a boundary term. The latter, of course, does 
not modify the local structure of quantum paths, hence the conclusion. 

A viable construction of quantum mechanics on fractional spacetimes opens up the 
possibility to test the latter as models of Nature. As stated in the introduction, eventu- 
ally one would like to make predictions in a field-theory context where the dimension of 
spacetime changes with the scale. However, at large-enough space/time scales the effect 
of dimensional flow is small and one can consider a fractional spacetime with fixed dimen- 
sionality close to 3+1, corresponding to an expansion in the small parameters e ~ 1 — ao,i- 
The scales at which this approximation holds are at least atomic [3], possibly smaller. 
It is in this context that the framework presented here can help to place bounds on the 

7 On the other hand, du cannot be greater than the topological dimension of the space wherein the curve 
is embedded. 
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theory from experiments in quantum mechanics. For example, the difference between the 
theoretical predictions and experimental data for the Lamb shift in hydrogen can be as- 
cribed to dimensional effects in order to give an upper bound for e. For theories with 
naive dimensional-regularization-like modifications, the constraint is |e| < 10 -11 at scales 
£ ~ 10 -11 m [56, 57]. It also turns out that this upper bound is more stringent than others 
at electroweak or astrophysical scales [3] . A revisit of the same problem in fractional space- 
times should be able to yield a similar constraint, probably of the same order of magnitude. 
This will be the subject of future investigations. 
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